Criterion method of GCCS for three-node VCSEL networks with delay coupling

ABSTRACT

A criterion method of GCCS (Globally Complete Chaos Synchronization) for three-node VCSEL (Vertical Cavity Surface Emitting Laser) networks with delay coupling is provided, including steps of: providing a delay-coupled VCSEL network consisting of three identical units and dynamic equations of the VCSEL network; providing assumptions of an outer-coupling matrix and a unitary matrix under the dynamic equations of the VCSEL network; in the three-node VCSEL network, determining rate equations of i-VCSEL, determining dynamic equations of a synchronization manifold, and determining a master-stability equation; calculating three maximum Lyapunov exponents; determining a stability of a synchronization state of the three-node VCSEL network, and determining whether the synchronization manifold of the VCSEL network is a chaotic waveform. Through a master-stability function, the method for determining whether the GCCS is achieved among all node lasers is provided, which solves a difficult problem of GCCS criterion for the VCSEL networks.

CROSS REFERENCE OF RELATED APPLICATION

The application claims priority under 35 U.S.C. 119(a-d) to CN 201810599588.5, filed Jun. 12, 2018.

BACKGROUND OF THE PRESENT INVENTION Field of Invention

The present invention relates to a field of VCSEL (Vertical Cavity Surface Emitting Laser) network, and more particularly to a criterion method of GCCS (Globally Complete Chaos Synchronization) for three-node VCSEL networks with delay coupling.

Description of Related Arts

The chaotic secure communication based on the semiconductor laser (SL) has distinct advantages, such as large bandwidth, high speed, low loss and high security. In future, it will have potential applications in the all-optical chaotic communication with the ultra-wide band and high speed. In the past twenty-four years, a series of achievements have been made in the development of chaotic secure communication system based on semiconductor lasers. Its development process for the topology is given as follows: point-to-point, point-to-multipoint, multipoint-to-multipoint, and complex network. Many researchers have explored the different topologies and the chaotic synchronization schemes of the SL-based secure communication. The researches about the topologies mainly focus on the delay connections between two nodes, and that among three or four nodes, which constitute the complex network topology. The GCCS (Globally Complete Chaos Synchronization) in these topologies plays an important role in complex network secure communication. However, the currently public GCCS only occur in a deterministic topology of three nodes or four nodes; and for three-node delay-coupled laser networks, there exist thirteen topology modules. If the semiconductor laser is subjected to optical feedback, there exist sixty-nine topology modules. For different topology modules of three-node lasers with delay coupling, the stability analysis of the synchronized manifold among three identical nodes has been a special focus due to the availability of analytic framework of the MSF (Master-Stability Function).

However, using the MSF, the criterion of the GCCS becomes complex if the node SL is subject to delay coupling, owing to the case that the delay-time could enhance and even induce synchronization in different topologies. Therefore, a lot of researches focus on the criterion of the GCCS in nonlinear system or complex network with delay coupling. For example, under large delay coupling, the calculation of the MSF for the synchronization criterion is obviously simplified; and under dual coupling channels, the complex behaviors of the GCCS are explored through utilizing the MSF. However, most of the previous researches about the MSF consider the condition that: the sum of each row in the outer-coupling matrix is zero (diffusive coupling), and the inner-coupling matrix is time-invariant. The researches about the GCCS of the complex network mainly consider the edge emitting laser as the node light source. However, in the complex chaotic secure communication network system, the VCSEL (Vertical Cavity Surface Emitting Laser) should serve as the node light source. Compared with the edge emitting laser, the VCSEL has many advantages, such as low threshold current, small volume, high efficiency, single-longitudinal-mode operation, and a circular output beam. However, the VCSEL also has more complex nonlinear dynamic behaviors and the unique polarization conversion characteristics. Besides the above conditions, the diversity of the topology modules of the three-node VCSEL network makes the criterion of the GCCS more complex. Therefore, the criterion of the GCCS in the VCSEL network is a key technical problem required to be solved urgently.

SUMMARY OF THE PRESENT INVENTION

An object of the present invention is to provide a criterion method of GCCS (Globally Complete Chaos Synchronization) for three-node VCSEL (Vertical Cavity Surface Emitting Laser) networks with delay coupling by using an MSF (Master-Stability Function) under nondiffusive and uniform delay coupling, so as to solve problems existing in the prior art.

The object of the present invention is realized through technical solutions as follows.

A criterion method of GCCS for three-node VCSEL networks with delay coupling comprises steps of:

(1), providing a delay-coupled VCSEL network consisting of three identical units, wherein dynamic equations of the VCSEL network are that:

${{\overset{g}{X_{i}}(t)} = {{F\left( {X_{i}(t)} \right)} + {\sum\limits_{j = 1}^{3}{K_{ij}{H\left( {{X_{i}(t)},{X_{j}\left( {t - \tau} \right)}} \right)}}}}},{i = 1},2,{3;}$

wherein: in the equations, X_(i)(t)=[X_(i1)(t), X_(i2)(t), X_(i3)(t), X_(i4)(t), X_(i5)(t), X_(i6)(t)]^(T)∈R⁶ represents a state variable of an i^(th) node; F:R⁶→R⁶ represents a nonlinear vector-valued function of an isolated node; H:R⁶×R⁶→R⁶ represents an inner-coupling function of the i^(th) node and a j^(th) node; τ is a coupling time-delay, which is assumed to be the same for all links; K=(K_(ij))₃₋₃∈R^(3×3) is an outer-coupling matrix which describes a coupling topology and a strength of each link in the network; K_(ii) is a self-feedback strength of the i^(th) node; and, K_(ij) is an injection strength from the j^(th) node to the i^(th) node;

(2), providing two assumptions, wherein:

for a first assumption, the sum of each row in the outer-coupling matrix K=(K_(ij))_(3×3)∈R^(3×3) is assumed to be a same constant (constant row sum for short) that:

${\sigma = {\sum\limits_{j = 1}^{3}K_{ij}}},{i = 1},2,{3;}$

for a second assumption, it is assumed that a unitary matrix U exists, which makes K^(T)=UΛU⁻¹;

the first assumption ensures an existence of an invariant synchronization manifold; and, for a given constant row sum σ, a dynamic in the synchronization manifold of the VCSEL network is that:

(t)=F(S(t))+σH(S(t),S(t−τ));

to a maximum Lyapunov exponent λ_(max) is a function of σ and λ_(k), and is calculated through a master-stability equation of:

(t)=Aξ(t)+λ_(k) DH(S(t−τ))ξ(t−τ),

in the equation, A=DF(S(t))+σDH(S(t)); λ_(k) is an eigenvalue of K, and k=1, 2, 3; DF and DH are Jacobian matrices calculated on the synchronization manifold; an MSF is able to be obtained through calculating the maximum Lyapunov exponent; the constant row sum σ is an eigenvalue of the outer-coupling matrix K; when λ_(k)=σ, the equation of

(t)=Aξ(t)+λ_(k)DH(S(t−τ))ξ(t−τ) is the master-stability equation for the synchronization manifold, and σ is related to disturbances within the synchronization manifold; and, transversal eigenvalues refer to all eigenvalues except for the eigenvalue σ;

(3), based on a spin-flip model and the dynamic equations of the VCSEL network, obtaining rate equations for an i-VCSEL in the three-node VCSEL networks that:

${\frac{{dE}_{x,y}^{i}(t)}{dt} = {{\left\{ {{\kappa\left\lbrack {{N^{i}(t)} - 1} \right\rbrack}m\;\gamma_{a}} \right\}{E_{x,y}^{i}(t)}} - {\kappa\;{n^{i}(t)}{E_{y,x}^{i}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack}}} \right\}} + {K_{i\; 1}{E_{x,y}^{1}\left( {t - \tau} \right)}{\cos\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{1}\left( {t - \tau} \right)}} \right\rbrack}} + {K_{i\; 2}{E_{x,y}^{2}\left( {t - \tau} \right)}{\cos\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{2}\left( {t - \tau} \right)}} \right\rbrack}} + {K_{i\; 3}{E_{x,y}^{3}\left( {t - \tau} \right)}{\cos\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{3}\left( {t - \tau} \right)}} \right\rbrack}}}},{\frac{d\;{\varphi_{x,y}^{i}(t)}}{dt} = {{{\kappa\;{\alpha\left\lbrack {{N^{i}(t)} - 1} \right\rbrack}m\;\gamma_{p}} \pm {\kappa\;{n^{i}(t)}\frac{E_{y,x}^{i}(t)}{E_{x,y}^{i}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack}} \right\}}} - {K_{i\; 1}\frac{E_{x,y}^{1}\left( {t - \tau} \right)}{E_{x,y}^{i}(t)}{\sin\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{1}\left( {t - \tau} \right)}} \right\rbrack}} - {K_{i\; 2}\frac{E_{x,y}^{2}\left( {t - \tau} \right)}{E_{x,y}^{i}(t)}{\sin\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{2}\left( {t - \tau} \right)}} \right\rbrack}} - {K_{i\; 3}\frac{E_{x,y}^{3}\left( {t - \tau} \right)}{E_{x,y}^{i}(t)}{\sin\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{3}\left( {t - \tau} \right)}} \right\rbrack}}}},{\frac{{dN}^{i}(t)}{dt} = {\gamma_{e}\left\{ {\mu - {{N^{i}(t)}\left\lbrack {1 + \left( {E_{x}^{i}(t)} \right)^{2} + \left( {E_{y}^{i}(t)} \right)^{2}} \right\rbrack} + {2\;{n^{i}(t)}{E_{x}^{i}(t)}{E_{y}^{i}(t)}{\sin\left( {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right)}}} \right\}}},{\frac{{dn}^{i}(t)}{dt} = {{{- \gamma_{s}}{n^{i}(t)}} - {\gamma_{e\;}\left\{ {{{n^{i}(t)}\left\lbrack {\left( {E_{x}^{i}(t)} \right)^{2} + \left( {E_{y}^{i}(t)} \right)^{2}} \right\rbrack} - {2\;{N^{i}(t)}{E_{x}^{i}(t)}{E_{y}^{i}(t)}{\sin\left( {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack}}} \right\}}}},$

wherein: in the equations, a superscript i represents the i-VCSEL; subscripts x and y respectively represent an x linear polarization mode and a y linear polarization mode; E is a slowly varying real amplitude of field; φ is a real phase of the field; N is a total carrier inversion between a conduction band and a valence band; n is a difference between carrier inversions for a spin-up radiation channel and a spin-down radiation channel; κ is a decay rate of the field; α is a linewidth enhancement factor; γ_(e) is a decay rate of total carrier population; γ_(s) is a spin-flip rate; γ_(a) and γ_(p) are linear anisotropies, respectively representing dichroism and birefringence; μ is a normalized injection current; and, central angular frequencies ω of all node VCSELs are assumed to be the same;

(4), according to the equation of

(t)=F(S(t))+σH(S(t),S(t−τ)), obtaining dynamic equations of the synchronization manifold of the three-node VCSEL network that:

${\frac{{dE}_{{sx},y}(t)}{dt} = {{\left\{ {{\kappa\left\lbrack {{N_{s}(t)} - 1} \right\rbrack}m\;\gamma_{a}} \right\}{E_{{sx},y}(t)}} - {\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}} + {\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}{\cos\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}}}},{\frac{d\;{\varphi_{{sx},y}(t)}}{dt} = {{\left\{ {\kappa\;{\alpha\left\lbrack {{N_{s}(t)} - 1} \right\rbrack}m\;\gamma_{p}} \right\} \pm {\frac{\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}}{E_{{sx},y}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}}} - {\frac{\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}}{E_{{sx},y}(t)}{\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}}}},{\frac{{dN}_{s}(t)}{dt} = {\gamma_{e}\left\{ {\mu - {{N_{s}(t)}\left\lbrack {1 + {E_{sx}^{2}(t)} + {E_{sy}^{2}(t)}} \right\rbrack} + {2\;{n_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}}},{\frac{{dn}_{s}(t)}{dt} = {{{- \gamma_{s}}{n_{s}(t)}} - {\gamma_{e\;}\left\{ {{{n_{s}(t)}\left\lbrack {{E_{sx}^{2}(t)} + {E_{sy}^{2}(t)}} \right\rbrack} - {2\;{N_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}}}},$

(5), according to the equation of

(t)=Aξ(t)+λ_(k)DH(S(t−τ))ξ(t−τ), obtaining master-stability equations for the three-node VCSEL network that:

${\frac{d\;\Delta\;{E_{x,y}(t)}}{dt} = {{\left\lbrack {{\kappa\left( {{N_{s}(t)} - 1} \right)}m\;\gamma_{a}} \right\rbrack\Delta\;{E_{x,y}(t)}} - {\kappa\;{n_{s}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{y,x}(t)}} - {\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}\left\{ {{- {\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \pm {\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{\varphi_{x}(t)}} - {\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}{\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{\varphi_{x,y}(t)}} - {\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{\varphi_{y}(t)}} + {\kappa\;{E_{{sx},y}(t)}\Delta\;{N(t)}} - {\kappa\;{E_{{sy},x}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{n(t)}} + {\lambda_{k}\begin{Bmatrix} {{{\cos\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{E_{x,y}\left( {t - \tau} \right)}} +} \\ {{E_{{sx},y}\left( {t - \tau} \right)}{\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{\varphi_{x,y}\left( {t - \tau} \right)}} \end{Bmatrix}}}},{\frac{d\;\Delta\;{\varphi_{x,y}(t)}}{dt} = {{\frac{m\;\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}}{E_{{sx},y}^{2}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{E_{x,y}(t)}} + {{{\frac{\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}}{E_{{sx},y}^{2}(t)}{\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{E_{x,y}(t)}} \pm {\frac{\kappa\;{n_{s}(t)}}{E_{{sx},y}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{E_{y,x}(t)}}} \pm {\frac{\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}}{E_{{sx},y}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{\varphi_{x}(t)}}} - {{\frac{\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}}{E_{{sx},y}(t)}{\cos\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right)}\Delta\;{\varphi_{x,y}(t)}} \pm {\frac{\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}}{E_{{sx},y}(t)}\left\{ {{- {\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}m\;\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{\varphi_{y}(t)}}} + {{\kappa\;\alpha\;\Delta\;{N(t)}} \pm {\frac{\kappa\;{E_{{sy},x}(t)}}{E_{{sx},y}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{n(t)}}} + {\lambda_{k}\left\{ {{\frac{- {\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}}{E_{{sx},y}(t)}\Delta\;{E_{x,y}\left( {t - \tau} \right)}} + {\frac{E_{{sx},y}\left( {t - \tau} \right)}{E_{{sx},y}(t)}{\cos\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{\varphi_{x,y}\left( {t - \tau} \right)}}} \right\}}}},{\frac{d\;\Delta\;{N(t)}}{dt} = {{\left\{ {{{- 2}\;\gamma_{e}{N_{s}(t)}{E_{sx}(t)}} + {2\;\gamma_{e}{n_{s}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{x}(t)}} + {\gamma_{e}\left\{ {{{- 2}\;{N_{s}(t)}{E_{sy}(t)}} + {2\;{n_{s}(t)}{E_{sx}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{y}(t)}} - {2\;\gamma_{e}{n_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{\varphi_{x}(t)}} + {2\;\gamma_{e}{n_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{\varphi_{y}(t)}} - {{\gamma_{e}\left\lbrack {1 + {E_{sx}^{2}(t)} + {E_{sy}^{2}(t)}} \right\rbrack}\Delta\;{N(t)}} + {2\;\gamma_{e}{E_{sx}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{n(t)}}}},{\frac{d\;\Delta\;{n(t)}}{dt} = {{\left\{ {{{- 2}\;\gamma_{e}{N_{s}(t)}{E_{sx}(t)}} + {2\;\gamma_{e}{n_{s}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{x}(t)}} + {\gamma_{e}\left\{ {{{- 2}\;{N_{s}(t)}{E_{sy}(t)}} + {2\;{n_{s}(t)}{E_{sx}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{y}(t)}} - {2\;\gamma_{e}{N_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{\varphi_{x}(t)}} + {2\;\gamma_{e}{N_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{\varphi_{y}(t)}} + {2\;\gamma_{e}{E_{sx}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{N(t)}} - {\left\{ {\gamma_{s} + {\gamma_{e}\left\lbrack {{E_{sx}^{2}(t)} + {E_{sy}^{2}(t)}} \right\rbrack}} \right\}\Delta\;{{n(t)}.}}}}$

(6), according to the equations in the steps (4) and (5), calculating three maximum Lyapunov exponents, determining a stability of a synchronization state of the three-node VCSEL network, and determining whether the synchronization manifold of the three-node VCSEL network is a chaotic waveform.

Preferably, in the step (1), when t→∞, if x₁(t)=x₂(t)=x₃(t)=S(t), the VCSEL network realizes globally complete synchronization, wherein S(t) represents the synchronization manifold of the VCSEL network; when S(t) is a chaotic waveform, the VCSEL network realizes the GCCS.

Preferably, in the step (2), for a given coupling topology with the constant row sum σ, if λ_(1 max) and λ_(2 max) are respectively negative at the transversal eigenvalues λ₁ and λ₂, the synchronization state is stable; if λ_(3 max) is positive at λ₃=σ, temporal traces of the synchronization manifold is chaotic.

The present invention has following beneficial effects.

For diverse three-node VCSEL networks with uniform time-delayed coupling, through the MSF, the present invention provides the method for determining whether the GCCS is achieved among all node lasers. The MSF can be obtained through calculating the MLE (Maximum Lyapunov Exponent) from the modified master-stability equation, and the MLE is a function of the constant row sum and the eigenvalue of the outer-coupling matrix. In the network, the outer-coupling matrix has the non-zero row sum (nondiffusive coupling), and the inner-coupling function is a function of self-node and delayed connection node. When two points determined by the constant row sum and two transversal eigenvalues are in the stability region (MLE is negative), the MLE as the function of the constant row sum and the eigenvalue (related to disturbances within the synchronization manifold) is positive, so that the arbitrarily given three-node VCSEL network is able to reach the GCCS, which solves the difficult problem of GCCS criterion in the VCSEL network.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a flow chart of a criterion method of GCCS (Globally Complete Chaos Synchronization) for three-node VCSEL (Vertical Cavity Surface Emitting Laser) networks with delay coupling according to the present invention.

FIG. 2 is a sketch view of an MSF (Master-Stability Function) for the three-node VCSEL network in a parameter space of σ and λ_(k) according to the present invention.

FIG. 3 is a sketch view of a ring topology based on three mutually coupled VCSELs according to the present invention.

FIG. 4 is a sketch view of detailed optical paths of the ring topology based on the three mutually coupled VCSELs according to the present invention.

FIG. 5 is a sketch view of temporal evolutions of synchronization errors for the three-node VCSEL network according to the present invention.

FIG. 6 is a sketch view of temporal traces of a synchronization manifold for the three-node VCSEL network according to the present invention.

FIG. 7 is a sketch view of a topology of a bus-shaped three-node VCSEL network according to the present invention.

FIG. 8 is a sketch view of detailed optical paths of the bus-shaped three-node VCSEL network according to the present invention.

FIG. 9 is a sketch view of temporal evolutions of synchronization errors for the bus-shaped three-node VCSEL network according to the present invention.

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT

Technical solutions of the present invention are further described in detail with accompanying drawings as follows, but the protection scope of the present invention is not limited thereto.

As shown in FIG. 1, a criterion method of GCCS (Globally Complete Chaos Synchronization) for three-node VCSEL (Vertical Cavity Surface Emitting Laser) networks with delay coupling comprises steps of:

(1), providing a delay-coupled VCSEL network consisting of three identical units, wherein dynamic equations of the VCSEL network are that:

${{{\overset{g}{X}}_{i}(t)} = {{F\left( {X_{i}(t)} \right)} + {\sum\limits_{j = 1}^{3}{K_{ij}{H\left( {{X_{i}(t)},{X_{j}\left( {t - \tau} \right)}} \right)}}}}},{i = 1},2,{3;}$

wherein: in the equations, X_(i)(t)=[X_(i1)(t), X_(i2)(t), X_(i3)(t), X_(i4)(t), X_(i5)(t), X_(i6)(t)]^(T)∈R⁶ represents a state variable of an i^(th) node; F:R⁶→R⁶ represents a nonlinear vector-valued function of an isolated node; H:R⁶×R⁶→R⁶ represents an inner-coupling function of the i^(th) node and a j^(th) node; r is a coupling time-delay, which is assumed to be the same for all links; K=(K_(ij))_(3×3)∈R^(3×3) is an outer-coupling matrix which describes a coupling topology and a strength of each link in the network; K_(ii) is a self-feedback strength of the i^(th) node; and, K_(ij) is an injection strength from the j^(th) node to the i^(th) node;

(2), providing two assumptions, wherein:

for a first assumption, the sum of each row in the outer-coupling matrix K=(K_(ij))_(3×3)∈R^(3×3) is assumed to be a same constant (constant row sum for short) that:

${\sigma = {\sum\limits_{j = 1}^{3}K_{ij}}},{i = 1},2,{3;}$

for a second assumption, it is assumed that a unitary matrix U exists, which makes K^(T)=UΛU⁻¹;

the first assumption ensures an existence of an invariant synchronization manifold; and, for a given constant row sum σ, a dynamic in the synchronization manifold of the VCSEL network is that:

(t)=F(S(t))+σH(S(t),S(t−τ)),

a maximum Lyapunov exponent λ_(max) is a function of σ and λ_(k), and is calculated through a master-stability equation of:

(t)=Aξ(t)+λ_(k) DH(S(t−τ))ξt(t−τ);

in the equation, A=DF(S(t))+σDH(S(t); λ_(k) is an eigenvalue of K, and k=1, 2, 3; DF and DH are Jacobian matrices calculated on the synchronization manifold; an MSF is able to be obtained through calculating the maximum Lyapunov exponent; the constant row sum σ is an eigenvalue of the outer-coupling matrix K; when λ_(k)=σ, the equation of

(t)=Aξ(t)+λ_(k)DH(S(t−τ))ξ(t−τ) is the master-stability equation for the synchronization manifold, and σ is related to disturbances within the synchronization manifold; and, transversal eigenvalues refer to all eigenvalues except for the eigenvalue σ;

(3), based on a spin-flip model and the dynamic equations of the VCSEL network, obtaining rate equations for an i-VCSEL in the three-node VCSEL networks that:

${\frac{{dE}_{x,y}^{i}(t)}{dt} = {{\left\{ {{\kappa\left\lbrack {{N^{i}(t)} - 1} \right\rbrack}m\;\gamma_{a}} \right\}{E_{x,y}^{i}(t)}} - {\kappa\;{n^{i}(t)}{E_{y,x}^{i}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack}}} \right\}} + {K_{i\; 1}{E_{x,y}^{1}\left( {t - \tau} \right)}{\cos\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{1}\left( {t - \tau} \right)}} \right\rbrack}} + {K_{i\; 2}{E_{x,y}^{2}\left( {t - \tau} \right)}{\cos\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{2}\left( {t - \tau} \right)}} \right\rbrack}} + {K_{i\; 3}{E_{x,y}^{3}\left( {t - \tau} \right)}{\cos\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{3}\left( {t - \tau} \right)}} \right\rbrack}}}},{\frac{d\;{\varphi_{x,y}^{i}(t)}}{dt} = {{{\kappa\;{\alpha\left\lbrack {{N^{i}(t)} - 1} \right\rbrack}m\;\gamma_{p}} \pm {\kappa\;{n^{i}(t)}\frac{E_{y,x}^{i}(t)}{E_{x,y}^{i}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack}} \right\}}} - {K_{i\; 1}\frac{E_{x,y}^{1}\left( {t - \tau} \right)}{E_{x,y}^{i}(t)}{\sin\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{1}\left( {t - \tau} \right)}} \right\rbrack}} - {K_{i\; 2}\frac{E_{x,y}^{2}\left( {t - \tau} \right)}{E_{x,y}^{i}(t)}{\sin\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{2}\left( {t - \tau} \right)}} \right\rbrack}} - {K_{i\; 3}\frac{E_{x,y}^{3}\left( {t - \tau} \right)}{E_{x,y}^{i}(t)}{\sin\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{3}\left( {t - \tau} \right)}} \right\rbrack}}}},{\frac{{dN}^{i}(t)}{dt} = {\gamma_{e}\left\{ {\mu - {{N^{i}(t)}\left\lbrack {1 + \left( {E_{x}^{i}(t)} \right)^{2} + \left( {E_{y}^{i}(t)} \right)^{2}} \right\rbrack} + {2\;{n^{i}(t)}{E_{x}^{i}(t)}{E_{y}^{i}(t)}{\sin\left( {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right)}}} \right\}}},{\frac{{dn}^{i}(t)}{dt} = {{{- \gamma_{s}}{n^{i}(t)}} - {\gamma_{e\;}\left\{ {{{n^{i}(t)}\left\lbrack {\left( {E_{x}^{i}(t)} \right)^{2} + \left( {E_{y}^{i}(t)} \right)^{2}} \right\rbrack} - {2\;{N^{i}(t)}{E_{x}^{i}(t)}{E_{y}^{i}(t)}{\sin\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack}}} \right\}}}},$

wherein: in the equations, a superscript i represents the i-VCSEL; subscripts x and y respectively represent an x linear polarization mode and a y linear polarization mode; E is a slowly varying real amplitude of field; φ is a real phase of the field; N is a total carrier inversion between a conduction band and a valence band; n is a difference between carrier inversions for a spin-up radiation channel and a spin-down radiation channel; κ is a decay rate of the field; α is a linewidth enhancement factor; γ_(e) is a decay rate of total carrier population; γ_(s) is a spin-flip rate; γ_(a) and γ_(p) are linear anisotropies, respectively representing dichroism and birefringence; μ is a normalized injection current; and, central angular frequencies ω of all node VCSELs are assumed to be the same;

(4), according to the equation of

(t)=F(S(t))+σH(S(t),S(t−τ)), obtaining dynamic equations of the synchronization manifold of the three-node VCSEL network that:

${\frac{{dE}_{{sx},y}(t)}{dt} = {{\left\{ {{\kappa\left\lbrack {{N_{s}(t)} - 1} \right\rbrack}m\;\gamma_{a}} \right\}{E_{{sx},y}(t)}} - {\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}} + {\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}{\cos\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}}}},{\frac{d\;{\varphi_{{sx},y}(t)}}{dt} = {{\left\{ {\kappa\;{\alpha\left\lbrack {{N_{s}(t)} - 1} \right\rbrack}m\;\gamma_{p}} \right\} \pm {\frac{\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}}{E_{{sx},y}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}}} - {\frac{\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}}{E_{{sx},y}(t)}{\sin\left\lbrack {{\sigma_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}}}},{\frac{{dN}_{s}(t)}{dt} = {\gamma_{e}\left\{ {\mu - {{N_{s}(t)}\left\lbrack {1 + {E_{sx}^{2}(t)} + {E_{sy}^{2}(t)}} \right\rbrack} + {2\;{n_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}}},{\frac{{dn}_{s}(t)}{dt} = {{{- \gamma_{s}}{n_{s}(t)}} - {\gamma_{e\;}\left\{ {{{n_{s}(t)}\left\lbrack {{E_{sx}^{2}(t)} + {E_{sx}^{2}(t)}} \right\rbrack} - {2\;{N_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}}}},$

(5), according to the equation of

(t)=Aξ(t)−λ_(k)DH(S(t−τ))ξ(t−τ), obtaining master-stability equations for the three-node VCSEL network that:

${\frac{d\;\Delta\;{E_{x,y}(t)}}{dt} = {{\left\lbrack {{\kappa\left( {{N_{s}(t)} - 1} \right)}m\;\gamma_{a}} \right\rbrack\Delta\;{E_{x,y}(t)}} - {\kappa\;{n_{s}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{y,x}(t)}} - {\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}\left\{ {{- {\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \pm {\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{\varphi_{x}(t)}} - {\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}{\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{\varphi_{x,y}(t)}} - {\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{\varphi_{y}(t)}} + {\kappa\;{E_{{sx},y}(t)}\Delta\;{N(t)}} - {\kappa\;{E_{{sy},x}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{n(t)}} + {\lambda_{k}\begin{Bmatrix} {{{\cos\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{E_{x,y}\left( {t - \tau} \right)}} +} \\ {{E_{{sx},y}\left( {t - \tau} \right)}{\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{\varphi_{x,y}\left( {t - \tau} \right)}} \end{Bmatrix}}}},{\frac{d\;\Delta\;{\varphi_{x,y}(t)}}{dt} = {{\frac{m\;\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}}{E_{{sx},y}^{2}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{E_{x,y}(t)}} + {{{\frac{\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}}{E_{{sx},y}^{2}(t)}{\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{E_{x,y}(t)}} \pm {\frac{\kappa\;{n_{s}(t)}}{E_{{sx},y}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{E_{y,x}(t)}}} \pm {\frac{\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}}{E_{{sx},y}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{\varphi_{x}(t)}}} - {{\frac{\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}}{E_{{sx},y}(t)}{\cos\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right)}\Delta\;{\varphi_{x,y}(t)}} \pm {\frac{\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}}{E_{{sx},y}(t)}\left\{ {{- {\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}m\;\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{\varphi_{y}(t)}}} + {{\kappa\;\alpha\;\Delta\;{N(t)}} \pm {\frac{\kappa\;{E_{{sy},x}(t)}}{E_{{sx},y}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{n(t)}}} + {\lambda_{k}\left\{ {{\frac{- {\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}}{E_{{sx},y}(t)}\Delta\;{E_{x,y}\left( {t - \tau} \right)}} + {\frac{E_{{sx},y}\left( {t - \tau} \right)}{E_{{sx},y}(t)}{\cos\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{\varphi_{x,y}\left( {t - \tau} \right)}}} \right\}}}},{\frac{d\;\Delta\;{N(t)}}{dt} = {{\left\{ {{{- 2}\;\gamma_{e}{N_{s}(t)}{E_{sx}(t)}} + {2\;\gamma_{e}{n_{s}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{x}(t)}} + {\gamma_{e}\left\{ {{{- 2}\;{N_{s}(t)}{E_{sy}(t)}} + {2\;{n_{s}(t)}{E_{sx}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{y}(t)}} - {2\;\gamma_{e}{n_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{\varphi_{x}(t)}} + {2\;\gamma_{e}{n_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{\varphi_{y}(t)}} - {{\gamma_{e}\left\lbrack {1 + {E_{sx}^{2}(t)} + {E_{sy}^{2}(t)}} \right\rbrack}\Delta\;{N(t)}} + {2\;\gamma_{e}{E_{sx}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{n(t)}}}},{\frac{d\;\Delta\;{n(t)}}{dt} = {{\left\{ {{{- 2}\;\gamma_{e}{N_{s}(t)}{E_{sx}(t)}} + {2\;\gamma_{e}{n_{s}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{x}(t)}} + {\gamma_{e}\left\{ {{{- 2}\;{N_{s}(t)}{E_{sy}(t)}} + {2\;{n_{s}(t)}{E_{sx}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{y}(t)}} - {2\;\gamma_{e}{N_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{\varphi_{x}(t)}} + {2\;\gamma_{e}{N_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{\varphi_{y}(t)}} + {2\;\gamma_{e}{E_{sx}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{N(t)}} - {\left\{ {\gamma_{s} + {\gamma_{e}\left\lbrack {{E_{sx}^{2}(t)} + {E_{sy}^{2}(t)}} \right\rbrack}} \right\}\Delta\;{{n(t)}.}}}}$

(6), according to the equations in the steps (4) and (5), calculating three maximum Lyapunov exponents, determining a stability of a synchronization state of the three-node VCSEL network, and determining whether the synchronization manifold of the three-node VCSEL network is a chaotic waveform.

Furthermore, in the step (1), when t→∞, if x₁(t)=x₂(t)=x₃(t)=S(t), the VCSEL network realizes globally complete synchronization, wherein S(t) represents the synchronization manifold of the VCSEL network; when S(t) is a chaotic waveform, the VCSEL network realizes the GCCS.

Furthermore, in the step (2), for a given coupling topology with the constant row sum σ, if λ_(1 max) and λ_(2 max) are respectively negative at the transversal eigenvalues λ₁ and λ₂, the synchronization state is stable; if λ_(3 max) is positive at λ₃=σ, temporal traces of the synchronization manifold is chaotic.

For three-node delay-coupled laser networks without the optical feedback, there exist thirteen possible topology modules, wherein: six out of the thirteen different modules are linear chains, while the other seven ones are ring structures. If the node VCSELs are subjected to an optical feedback injection, the possible module number is sixty-nine. In order to explore the GCCS properties for different topologies in the three-node VCSEL network, the MLE (Maximum Lyapunov Exponent) is required to be calculated. The MLE is a function of the constant row sum σ and the eigenvalue λ_(k), and parameter values for numerical simulation are showed as follows.

Parameter Symbol Value Linewidth enhancement factor α 3 Dichroism γ_(a)/ns⁻¹ −0.1 Decay rate of carrier population γ_(e)/ns⁻¹ 1 Normalized injection current μ 1.5 Decay rate of field κ/ns⁻¹ 300 Birefringence γ_(p)/ns⁻¹ 2 Spin-flip rate γ_(s)/ns⁻¹ 50 Central wavelength λ/nm 850 Time-delay τ/ns 4

For simplicity, the real eigenvalue case is mainly considered herein. Specifically, the constant row sum σ∈[10¹⁰, 6×10¹⁰] and the eigenvalue λ_(k)∈[−3×10¹⁰, 6×10¹⁰] are considered; and FIG. 2 shows the MSF for the three-node VCSEL network in a parameter space of σ and λ_(k). It can be seen from FIG. 2 that: a stability region is determined by λ_(max)<0; and, when λ_(k)=σ, the corresponding λ_(max) is positive, indicating that all of the synchronization manifolds are chaotic waveforms.

Then, two given topology modules are respectively taken as two examples, so as to explore the validity of the MSF in FIG. 2. It should be noted that: the two fixed topologies are required to meet the first and second assumptions, and initial conditions of three VCSELs are slightly different. In order to describe the stability of the synchronization state, synchronization errors among three node-VCSELs are defined that: e ₁₂(t)=X ₁(t)−X ₂(t)=[e ₁₂₁(t),e ₁₂₂(t),e ₁₂₃(t),e ₁₂₄(t),e ₁₂₅(t),e ₁₂₆(t)]^(T) e ₂₃(t)=X ₂(t)−X ₃(t)=[e ₂₃₁(t),e ₂₃₂(t),e ₂₃₃(t),e ₂₃₄(t),e ₂₃₅(t),e ₂₃₆(t)]^(T) e ₁₃(t)=X ₁(t)−X ₃(t)=[e ₁₃₁(t),e ₁₃₂(t),e ₁₃₃(t),e ₁₃₄(t),e ₁₃₅(t),e ₁₃₆(t)]_(T);

wherein: X_(i)(t) are the state variables in the step (1); for e_(ijq)(t), i and j respectively represent i-VCSEL and j-VCSEL, q is a dimension of the state variables of the node-VCSEL, and q=1, 2, 3, 4, 5, 6 respectively represent E_(x), E_(y), φ_(x), φ_(y), N and n.

The ring topology with three mutually coupled node VCSELs is showed in FIG. 3. FIG. 4 shows detailed optical paths for performing the GCCS, wherein the OIs (Optical Isolators), respectively OI₁, OI₂ and OI₃, are for guaranteeing the unidirectional propagation of the optical wave. Mutual coupling among three node-VCSELs is realized through BSs (Beam Splitters), BS₁-BS₅. The node-VCSEL1 and the node-VCSEL2 are subjected to the optical feedback respectively through the M₂ (Mirror 2) and the M₅ (Mirror 5). When the bias current, the optical feedback strength and the injection strength are fixed to certain values, all of the node-VCSELs emit two chaotic PCs (Polarization Components). Because polarization directions thereof are along directions of an x-axis and a y-axis, the PCs are respectively defined as X-LP (X-Linear Polarization) and Y-LP (Y-Linear Polarization). For the i(j)-node VCSEL, the two PCs thereof are respectively denoted by X_(i(j))-LP and Y_(i(j))-LP, and decomposed by the PBS (Polarization Beam Splitter) with the number of i(j); after being converted into the electric signals by the PD (Photodetector), the temporal waveforms of the X_(i)-LP (Y_(i)-LP) and X_(j)-LP (Y_(j)-LP) may be synchronized to each other, under i≠j. The synchronization error e_(ij1)(t) is able to be obtained through the synchronous subtraction between the X_(i)-LP and the X_(j)-LP. Similarly, e_(ij2)(t) is able to be obtained through the synchronous subtraction between the Y_(i)-LP and the Y_(j)-LP. The outer-coupling matrix herein is set as:

$K = {\begin{bmatrix} 10^{10} & 10^{10} & 10^{10} \\ 10^{10} & 10^{10} & 10^{10} \\ {2 \times 10^{10}} & 10^{10} & 0 \end{bmatrix}.}$

A set of eigenvalues of K are that: λ₁=1×10¹⁰, λ₂=−6.2867×10⁻⁷, and λ₃=σ=3×10¹⁰. When σ=3×10¹⁰, it can be seen from FIG. 2 that: in the parameter space of σ and λ_(k); the points (σ, λ₁) and (σ, λ₂) are both in the stability region, indicating that the complete synchronization can be achieved among three node-VCSELs; λ_(max) corresponding to the point (σ, λ₃) is positive, indicating that the synchronization manifold is a chaotic waveform. In order to verify the validity of the above results, FIG. 5 shows temporal evolutions of the synchronization errors for the three-node VCSEL network. It can be seen from FIG. 5 that: although the initial value of each VCSEL is different, all of the error curves violently oscillate at the initial stage of the evolution, then tend to zero rapidly, and finally keep at zero. All of the curves tend to zero after the transient state (about 100 ns), indicating that the globally asymptotic synchronization among three node-VCSELs can be achieved. After reaching the complete synchronization, FIG. 6 further shows the temporal traces of the synchronization manifold for the three-node VCSEL network. It can be obviously seen from FIG. 6 that: the synchronization manifold exhibits the chaotic behaviors. Therefore, the above results completely verify the GCCS criterion by the MSF displayed in FIG. 2.

A bus-shaped three-node VCSEL network is taken as the other example to further verify the validity of the GCCS criterion. FIG. 7 is a sketch view of a topology of the bus-shaped three-node VCSEL network. FIG. 8 shows detailed optical paths for the bus-shaped three-node VCSEL network. In FIGS. 7 and 8, the VCSEL2 is subjected to an optical feedback through the M₂; the VCSEL3 is subject to the unidirectional light injection from the VCSEL2 output; mutual coupling between the VCSEL1 and the VCSEL2 is realized through the BS₁ and the BS₂; the synchronization errors e_(ij1)(t) and e_(ij2)(t) can be obtained through the same method in FIG. 4. For the above configuration, the outer-coupling matrix is set to be:

$K = {\begin{bmatrix} 0 & {2 \times 10^{10}} & 0 \\ 10^{10} & 10^{10} & 0 \\ 0 & {2 \times 10^{10}} & 0 \end{bmatrix}.}$

A set of eigenvalues of K are that: λ₁=−1×10¹⁰, λ₂=0, and λ₃=σ=2×10¹⁰. When σ=2×10¹⁰, it can be seen from FIG. 2 that: the point (σ, λ₁) is not in the stability region, and the point (σ, λ₂) is in the stability region, indicating that the topology is unable to reach the globally asymptotic synchronization. In order to further observe the synchronization properties among three-node lasers, FIG. 9 shows temporal evolutions of the synchronization errors for the bus-shaped three-node VCSEL network. It can be seen from FIG. 9 that: both of e_(12q)(t) and e_(23q)(t) violently oscillate during the whole evolution process, indicating that the synchronization between the VCSEL1 and the VCSEL2, as well as that between the VCSEL2 and the VCSEL3 cannot be achieved. Meanwhile, e_(13q)(t) have similar temporal evolutions to those given in FIG. 5, indicating that the complete synchronization can be achieved between the VCSEL1 and the VCSEL3. From above results, it can be concluded that the global synchronization among three-node lasers are impossible to be achieved in the bus-shaped topology network. The above-mentioned results follow closely the theoretical prediction, further indicating the effectiveness of the theory.

It should be illustrated that the above-described is merely the preferred embodiment of the present invention. It should be also noted that the various improvements and modifications made by one skilled in the art without departing from the principle of the present invention are all encompassed in the protection scope of the present invention. 

What is claimed is:
 1. A criterion method of GCCS (Globally Complete Chaos Synchronization) for three-node VCSEL (Vertical Cavity Surface Emitting Laser) networks with delay coupling, comprising steps of: (1), providing a delay-coupled VCSEL network consisting of three identical units, which comprises specific steps of: using multiple optical isolators in optical paths for performing the GCCS, so as to guarantee unidirectional propagation of optical wave; processing the three identical units with mutual coupling through multiple beam splitters; providing optical feedbacks to a first VCSEL node and a second VCSEL node through a first mirror and a second minor, respectively; wherein when a bias current, an optical feedback strength and an injection strength are fixed to certain values, all of the three identical units emit two chaotic polarization components; because polarization directions thereof are along directions of an x-axis and a y-axis, the chaotic polarization components are defined as x linear polarization and y linear polarization; dynamic equations of the VCSEL network are that: ${{{\overset{g}{X}}_{i}(t)} = {{F\left( {X_{i}(t)} \right)} + {\sum\limits_{j = 1}^{3}{K_{ij}{H\left( {{X_{i}(t)},{X_{j}\left( {t - \tau} \right)}} \right)}}}}},{i = 1},2,{3;}$ wherein: in the equations, X_(i) (t)=[X_(i1)(t), X_(i2) (t), X_(i3) (t), X_(i4) (t), X_(i5)(t), X_(i6) (t)]^(T) ∈R⁶ represents a state variable of an i^(th) node; F: R⁶→R⁶ represents a nonlinear vector-valued function of an isolated node; H: R⁶×R⁶→R⁶ represents an inner-coupling function of the i^(th) node and a j^(th) node; τ is a coupling time-delay, which is assumed to be the same for all links; K=(K_(ij))_(3×3) ∈R^(3×3) is an outer-coupling matrix which describes a coupling topology and a strength of each link in the network; K_(ii) is a self-feedback strength of the i^(th) node; and, K_(ij) is an injection strength from the j^(th) node to the i^(th) node; (2), providing two assumptions, wherein: for a first assumption, the sum of each row in the outer-coupling matrix K=(K_(ij))_(3×3) ∈R^(3×3) is assumed to be a same constant (constant row sum for short) that: ${\sigma = {\sum\limits_{j = 1}^{3}K_{ij}}},{i = 1},2,{3;}$ for a second assumption, it is assumed that a unitary matrix U exists, which makes K^(T)=UΛU⁻¹; the first assumption ensures an existence of an invariant synchronization manifold; and, for a given constant row sum σ, a dynamic in the synchronization manifold of the VCSEL network is that:

(t)=F(S(t))+σH(S(t),S(t−τ)); a maximum Lyapunov exponent λ_(max) is a function of σ and λ_(k), and is calculated through a master-stability equation of:

(t)=Aξ(t)+λ_(k) DH(S(t−τ))ξ(t−τ); in the equation, A=DF(S(t))+σDH (S(t)); λ_(k) is an eigenvalue of K, and k=1, 2, 3; DF and DH are Jacobian matrices calculated on the synchronization manifold; an MSF (Master-Stability Function) is able to be obtained through calculating the maximum Lyapunov exponent; the constant row sum σ is an eigenvalue of the outer-coupling matrix K; when λ_(k)=σ, the equation of

(t)=Aξ(t)+λ_(k)DH(S(t−τ))ξ(t−τ) is the master-stability equation for the synchronization manifold, and σ is related to disturbances within the synchronization manifold; and, transversal eigenvalues refer to all eigenvalues except for the eigenvalue σ; (3), based on a spin-flip model and the dynamic equations of the VCSEL network, obtaining rate equations for an i-VCSEL in the three-node VCSEL networks that: ${\frac{{dE}_{x,y}^{i}(t)}{dt} = {{\left\{ {{\kappa\left\lbrack {{N^{i}(t)} - 1} \right\rbrack}m\;\gamma_{a}} \right\}{E_{x.y}^{i}(t)}} - {\kappa\;{n^{i}(t)}{E_{y,x}^{i}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack}}} \right\}} + {K_{i\; 1}{E_{x,y}^{1}\left( {t - \tau} \right)}{\cos\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{1}\left( {t - \tau} \right)}} \right\rbrack}} + {K_{i\; 2}{E_{x,y}^{2}\left( {t - \tau} \right)}{\cos\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{2}\left( {t - \tau} \right)}} \right\rbrack}} + {K_{i\; 3}{E_{x,y}^{3}\left( {t - \tau} \right)}{\cos\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{3}\left( {t - \tau} \right)}} \right\rbrack}}}},{\frac{d\;{\varphi_{x,y}^{i}(t)}}{dt} = {{{\kappa\;{\alpha\left\lbrack {{N^{i}(t)} - 1} \right\rbrack}m\;\gamma_{p}} \pm {\kappa\;{n^{i}(t)}\frac{E_{y,x}^{i}(t)}{E_{x,y}^{i}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack}} \right\}}} - {K_{i\; 1}\frac{E_{x,y}^{1}\left( {t - \tau} \right)}{E_{x,y}^{i}(t)}{\sin\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{1}\left( {t - \tau} \right)}} \right\rbrack}} - {K_{i\; 2}\frac{E_{x,y}^{2}\left( {t - \tau} \right)}{E_{x,y}^{i}(t)}{\sin\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{2}\left( {t - \tau} \right)}} \right\rbrack}} - {K_{i\; 3}\frac{E_{x,y}^{3}\left( {t - \tau} \right)}{E_{x,y}^{i}(t)}{\sin\left\lbrack {{\omega\;\tau} + {\varphi_{x,y}^{i}(t)} - {\varphi_{x,y}^{3}\left( {t - \tau} \right)}} \right\rbrack}}}},{\frac{{dN}^{i}(t)}{dt} = {\gamma_{e}\left\{ {\mu - {{N^{i}(t)}\left\lbrack {1 + \left( {E_{x}^{i}(t)} \right)^{2} + \left( {E_{y}^{i}(t)} \right)^{2}} \right\rbrack} + {2\;{n^{i}(t)}{E_{x}^{i}(t)}{E_{y}^{i}(t)}{\sin\left( {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right)}}} \right\}}},{\frac{{dn}^{i}(t)}{dt} = {{{- \gamma_{s}}{n^{i}(t)}} - {\gamma_{e\;}\left\{ {{{n^{i}(t)}\left\lbrack {\left( {E_{x}^{i}(t)} \right)^{2} + \left( {E_{y}^{i}(t)} \right)^{2}} \right\rbrack} - {2\;{N^{i}(t)}{E_{x}^{i}(t)}{E_{y}^{i}(t)}{\sin\left\lbrack {{\varphi_{y}^{i}(t)} - {\varphi_{x}^{i}(t)}} \right\rbrack}}} \right\}}}},$ wherein: in the equations, a superscript i represents the i-VCSEL; subscripts x and y respectively represent an x linear polarization mode and a y linear polarization mode; E is a slowly varying real amplitude of field; φ is a real phase of the field; N is a total carrier inversion between a conduction band and a valence band; n is a difference between carrier inversions for a spin-up radiation channel and a spin-down radiation channel; κ is a decay rate of the field; α is a linewidth enhancement factor; γ_(e) is a decay rate of total carrier population; γ_(s) is a spin-flip rate; γ_(a) and γ_(p) are linear anisotropies, respectively representing dichroism and birefringence; μ is a normalized injection current; and, central angular frequencies ω of all node VCSELs are assumed to be the same; (4), according to the equation of

(t)=F(S(t))+σH(S(t), S(t−τ)), obtaining dynamic equations of the synchronization manifold of the three-node VCSEL network that: ${\frac{{dE}_{{sx},y}(t)}{dt} = {{\left\{ {{\kappa\left\lbrack {{N_{s}(t)} - 1} \right\rbrack}m\;\gamma_{a}} \right\}{E_{{sx},y}(t)}} - {\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}} + {\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}{\cos\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}}}},{\frac{d\;{\varphi_{{sx},y}(t)}}{dt} = {{\left\{ {\kappa\;{\alpha\left\lbrack {{N_{s}(t)} - 1} \right\rbrack}m\;\gamma_{p}} \right\} \pm {\frac{\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}}{E_{{sx},y}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}}} - {\frac{\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}}{E_{{sx},y}(t)}{\sin\left\lbrack {{\sigma_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}}}},{\frac{{dN}_{s}(t)}{dt} = {\gamma_{e}\left\{ {\mu - {{N_{s}(t)}\left\lbrack {1 + {E_{sx}^{2}(t)} + {E_{sy}^{2}(t)}} \right\rbrack} + {2\;{n_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}}},{\frac{{dn}_{s}(t)}{dt} = {{{- \gamma_{s}}{n_{s}(t)}} - {\gamma_{e\;}\left\{ {{{n_{s}(t)}\left\lbrack {{E_{sx}^{2}(t)} + {E_{sy}^{2}(t)}} \right\rbrack} - {2\;{N_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}}}},$ (5), according to the equation of

(t)=Aξ(t)+λ_(k)DH(S(t−τ))ξ(t−τ), obtaining master-stability equations for the three-node VCSEL network that: ${\frac{d\;\Delta\;{E_{x,y}(t)}}{dt} = {{\left\lbrack {{\kappa\left( {{N_{s}(t)} - 1} \right)}m\;\gamma_{a}} \right\rbrack\Delta\;{E_{x,y}(t)}} - {\kappa\;{n_{s}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{y,x}(t)}} - {\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}\left\{ {{- {\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \pm {\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{\varphi_{x}(t)}} - {\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}{\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{\varphi_{x,y}(t)}} - {\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{\varphi_{y}(t)}} + {\kappa\;{E_{{sx},y}(t)}\Delta\;{N(t)}} - {\kappa\;{E_{{sy},x}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{n(t)}} + {\lambda_{k}\begin{Bmatrix} {{{\cos\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{E_{x,y}\left( {t - \tau} \right)}} +} \\ {{E_{{sx},y}\left( {t - \tau} \right)}{\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{\varphi_{x,y}\left( {t - \tau} \right)}} \end{Bmatrix}}}},{\frac{d\;\Delta\;{\varphi_{x,y}(t)}}{dt} = {{\frac{m\;\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}}{E_{{sx},y}^{2}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{E_{x,y}(t)}} + {{{\frac{\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}}{E_{{sx},y}^{2}(t)}{\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{E_{x,y}(t)}} \pm {\frac{\kappa\;{n_{s}(t)}}{E_{{sx},y}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{E_{y,x}(t)}}} \pm {\frac{\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}}{E_{{sx},y}(t)}\left\{ {{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack} \pm {\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{\varphi_{x}(t)}}} - {{\frac{\sigma\;{E_{{sx},y}\left( {t - \tau} \right)}}{E_{{sx},y}(t)}{\cos\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right)}\Delta\;{\varphi_{x,y}(t)}} \pm {\frac{\kappa\;{n_{s}(t)}{E_{{sy},x}(t)}}{E_{{sx},y}(t)}\left\{ {{- {\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}m\;\alpha\;{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{\varphi_{y}(t)}}} + {{\kappa\;\alpha\;\Delta\;{N(t)}} \pm {\frac{\kappa\;{E_{{sy},x}(t)}}{E_{{sx},y}(t)}\left\{ {{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}m\;\alpha\;{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}} \right\}\Delta\;{n(t)}}} + {\lambda_{k}\left\{ {{\frac{- {\sin\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}}{E_{{sx},y}(t)}\Delta\;{E_{x,y}\left( {t - \tau} \right)}} + {\frac{E_{{sx},y}\left( {t - \tau} \right)}{E_{{sx},y}(t)}{\cos\left\lbrack {{\varphi_{{sx},y}(t)} - {\varphi_{{sx},y}\left( {t - \tau} \right)} + {\omega\;\tau}} \right\rbrack}\Delta\;{\varphi_{x,y}\left( {t - \tau} \right)}}} \right\}}}},{\frac{d\;\Delta\;{N(t)}}{dt} = {{\left\{ {{{- 2}\;\gamma_{e}{N_{s}(t)}{E_{sx}(t)}} + {2\;\gamma_{e}{n_{s}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{x}(t)}} + {\gamma_{e}\left\{ {{{- 2}\;{N_{s}(t)}{E_{sy}(t)}} + {2\;{n_{s}(t)}{E_{sx}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{y}(t)}} - {2\;\gamma_{e}{n_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{\varphi_{x}(t)}} + {2\;\gamma_{e}{n_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{\varphi_{y}(t)}} - {{\gamma_{e}\left\lbrack {1 + {E_{sx}^{2}(t)} + {E_{sy}^{2}(t)}} \right\rbrack}\Delta\;{N(t)}} + {2\;\gamma_{e}{E_{sx}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{n(t)}}}},{\frac{d\;\Delta\;{n(t)}}{dt} = {{\left\{ {{{- 2}\;\gamma_{e}{N_{s}(t)}{E_{sx}(t)}} + {2\;\gamma_{e}{n_{s}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{x}(t)}} + {\gamma_{e}\left\{ {{{- 2}\;{N_{s}(t)}{E_{sy}(t)}} + {2\;{n_{s}(t)}{E_{sx}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}}} \right\}\Delta\;{E_{y}(t)}} - {2\;\gamma_{e}{N_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{\varphi_{x}(t)}} + {2\;\gamma_{e}{N_{s}(t)}{E_{sx}(t)}{E_{sy}(t)}{\cos\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{\varphi_{y}(t)}} + {2\;\gamma_{e}{E_{sx}(t)}{E_{sy}(t)}{\sin\left\lbrack {{\varphi_{sy}(t)} - {\varphi_{sx}(t)}} \right\rbrack}\Delta\;{N(t)}} - {\left\{ {\gamma_{s} + {\gamma_{e}\left\lbrack {{E_{sx}^{2}(t)} + {E_{sy}^{2}(t)}} \right\rbrack}} \right\}\Delta\;{{n(t)}.}}}}$ (6), according to the equations in the steps (4) and (5), calculating three maximum Lyapunov exponents, determining a stability of a synchronization state of the three-node VCSEL network, and determining whether the synchronization manifold of the three-node VCSEL network is a chaotic waveform.
 2. The criterion method of the GCCS for the three-node VCSEL networks with the delay coupling, as recited in claim 1, wherein: in the step (1), when t→∞, if x₁(t)=x₂(t)=x₃(t)=S(t), the VCSEL network realizes globally complete synchronization, wherein S(t) represents the synchronization manifold of the VCSEL network; when S(t) is a chaotic waveform, the VCSEL network realizes the GCCS.
 3. The criterion method of the GCCS for the three-node VCSEL networks with the delay coupling, as recited in claim 1, wherein: in the step (2), for a given coupling topology with the constant row sum σ, if λ_(1 max) and λ_(2 max) are respectively negative at the transversal eigenvalues λ₁ and λ₂ the synchronization state is stable; if λ_(3 max) is positive at λ₃=σ, temporal traces of the synchronization manifold is chaotic. 